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ABSTRACT 

The presence of unburned material in the ejecta of normal Type la supernovae (SNe la) is investigated 
using early-time spectroscopy obtained by the Carnegie Supernova Project (CSP). The tell-tale signa- 
ture of pristine material from a C+O white dwarf progenitor star is the presence of carbon, as oxygen 
is also a product of carbon burning. The most prominent carbon lines in optical spectra of SNe la are 
expected to arise from C II. We find that at least 30% of the objects in the sample show an absorption 
at «6300 A which is attributed to C II A6580. An alternative identification of this absorption as Ha 
is considered to be unlikely. These findings imply a larger incidence of carbon in SNe la ejecta than 
previously noted. We show how observational biases and physical conditions may hide the presence 
of weak C II lines, and account for the scarcity of previous carbon detections in the literature. This 
relatively large frequency of carbon detections has crucial implications on our understanding of the 
explosive process. Furthermore, the identification of the 6300 A absorptions as carbon would imply 
that unburned material is present at very low expansion velocities, merely wlOOO km s _1 above the 
bulk of Si II. Based on spectral modeling, it is found that the detections are consistent with a mass of 
carbon of 10 -3 - 10 -2 M Q . The presence of this material so deep in the ejecta would imply substantial 
mixing, which may be related to asymmetries of the flame propagation. Another possible explanation 
for the carbon absorptions may be the existence of clumps of unburned material along the line of sight. 
However, the uniformity of the relation between C II and Si II velocities is not consistent with such 
small-scale asymmetries. The spectroscopic and photometric properties of SNe la with and without 
carbon signatures are compared. A trend toward bluer color and lower luminosity at maximum light 
is found for objects which show carbon. 

Subject headings: supernovae: general - techniques: spectroscopic 
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1. INTRODUCTION 

The generally favored picture of a Type la SN (SN la) 
is the thermonuclear expl osion of a C+O white dwarf 
(WD) star as proposed by iHovle fc Fowler] ([1960) . The 
WD is assumed to belong to a binary system, and to 
explode when it reaches the Chandrasekhar mass limit 
(Mch) via two main, alternative channels. In the single- 
degenerate scenario, a main-sequence or giant companion 
star transfers mass to t he WD by Ro c he-Lobe overflow 
(IWhelan fc Ibenl H983I: INomotol ll982T: llben fc Tutukovl 
11984) . Alternatively, in the double-degenerate scenario, 
the companion star may be another WD, and the Mch 
limit is reached or exceeded when both objects coalesce 
(|Iben fc Tutukovl Il984t IWebbinkl 119841 ). The actual in- 
cidence of each of these evolutionary paths is currently 
under debate. 

Further controversy regarding the physical nature of 
SNe la is related to the mechanism by which the 
thermonuclear flame propaga tes inside the WD (see 
iHillebrandt fc Niemeverl 120001 and references therein) . 
It is generally believed the e xplosion initiates as a de- 
flagra tion, or subsonic flame (jNomoto. Sugimoto fc Neol 
1976). If the WD explodes as a prompt detonation, then 
no intermediate-mass elements would be left in the outer 
ejecta, as opposed to the observational evidence, and 
in contradiction w ith the chemical evolution of the in- 
terstellar medium (jArnett. Truran fc Wooslevl 119711 ). A 
pure deflagration, however, would not provide the nec- 
essary energy and n ucleosynthesis to explain most nor - 
mal SNe la (see, e.g. iGamezoT Khokhlov fc~Q7anl[2005h . 
It seems that a delayed detonation, i.e. a transition 
to a supersonic regime after the WD is pre-expanded 
in the deflagration phase, is necessary in order to re- 
produce the typical and bright SNe la (jKhokhlovl 119911 : 
IHoflich. Khokhlov fc Wheeler! I1995D . The mechanism 
that produces the transition is not completely under- 
stood, mostly due to the complexity of flame-propagation 
physics and the limited resolution of current models. 

Despite these unknowns, SNe la are successfully em- 
ployed as precision cosmological distance indicators due 
to their homogeneous properties. Their peak luminosi- 
ties can be calibrated through the observed light-curve 
width and color to provide di stances with a precisi on of 
^0.1 mag or better (e.g., see iFolatelli et al.ll2010bl and 
references therein). This calibrated precision led to the 
discovery of the acc eleration in the current expansion 
rate of the Universe (|Riess et al.lll998t iPerlmutter et al.l 
[19991) . 

Some spectral characteristics, such as the strength of 
Si II lines near maximum light, show variations that 
are correlated with photometric properties. Temperature 
in the ejecta has been identified as a pos sible driver of 
this l uminosity-related spectral diversity (jNugent et al.l 
1995). Temperature in turn depends on the amount 
of 56 Ni synthesized during the explosion, thus provid- 
ing a neat physical picture for SNe la. However, this 
one-parameter picture is insufficient to describe inhomo- 
geneities in other photometric and spectroscopic prop- 
erties that have been observed to vary independently 
of luminosity. Moreover, the connection of 56 Ni mass 
with temperature produces changes in the opacity which 
lead to more complex ef fects in the observed proper - 
ties (|Hoflich fc Khokhlov|[l99l iPinto fc Eastmanll2000f) . 



iKasen fc Wooslevl (|2007l ) performed a detailed study of 
the line blanketing effects on the light curves and spectra. 

The origin of these inhomogeneities has recently 
been linked to varying properties of the progenitor 
and explosion mec hanism (e.g., see IHoflich et alj 120101 : 
IMaeda et al.l 120111 ) . The latter determines the extent 
of burning and the degree of mixing of ashes and un- 
burned fuel within the ejecta. During the last decade, 
multi-dimensional calculations have been introduced in 
an attempt to fully reproduce the burni ng process from 
first principles (see IMaeda et al.l I2010bl and references 
therein). These models, although still relatively rudi- 
mentary, naturally introduce a variety of effects which 
may succeed in reproducing several aspects of the ob- 
served spectral diversity. For example, an asymmetric 
explosion, in which the thermonuclear reaction is ignited 
at an offset from the WD center, provides an explana- 
tion not only for the obser ved scatter of the peak mag- 
nitudes (|Kasen et al.l 120091) . but also for the diverse ve- 
locity and velocity g radient seen in the Si II A6355 line 
(|Maeda et al.ll2010d ). 

It is an observational fact that spectral inhomo- 
geneities are stronger before o r near maximum ligh t than 
later on ()Benetti et al.l I2005t IBranch et all 120071. At 
these early stages the outer part of the ejecta is ob- 
served. Under the as sumption of homologous expansion 
(e.g., see Rdpke 2005|), this outer material moves at high 
velocities (v >10000 km s _1 ). As a consequence, obser- 
vations indicate that, whatever the variations in the pro- 
genitor system and explosion mechanism are, the most 
noticeable effects appear in the fast-expanding, outer re- 
gions of the ejecta. Variations are commonly observed 
among normal SNe la in the relative strength of spectral 
lines, the evolution of expansion velocities, the appear- 
ance of high-velocity components, and the presence of 
dif ferent ionization states, among other properties (e.g., 
seelNugent et al.lll995HBenetti et al.H2005HMazzali et alj 
120051 ). The study of pre-maximum spectra thus promises 
to be very useful in the validation of proposed models. 

A relevant characteristic which can be studied using 
early-time spectra is the possible presence of unburned 
material in the outer ejecta. Most of the available mod- 
els predict that part of the primordial WD material 
composed of carbon and oxygen is left unburned. The 
amount and location of this unburned material has direct 
implications in the validity of the models. Since oxygen 
is also a product of carbon burning, it is the detection of 
carbon what would most clearly evidence the existence of 
unburned material. Furthermore, singly-ionized states of 
both C and O are expecte d to be dominant in the outer 
ejecta (v >10000 km s" 1 : iTanaka et al.1 120081 ). Due to 
high excitation temperatures of O II, however, no strong 
lines of this ion appear in the optical range. Some strong 
O I transitions produce observed absorptions which are 
blended with lines of Mg II, another product of carbon 
burning. At sufficiently low temperature, C I lines are 
expected to appear , most prominently in the near-IR 
(M arion et aLll2006D . however this condition is not likely 
to occur in the regions of the ejecta probed before max- 
imum light. The same is true for the high-temperature 
condition required to produce C III lines, unless non- 
LTE effects play an important role. Consequently, the 
possible detection of unburned material must rely on the 
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presence of C II lines. 

Until recently, the detection of carbon lines in SN-Ia 
spectra was uncommon. C II lines have been clearly 
identified in optical spectra of a few peculiarly luminous, 
slowly expanding objects which have been suggested to 
arise from super- Chandrasekhar mass progenitors (e.g., 
Howell et all 120061: jHicken et al.l [2007t lYamanaka et all 
20091* IScalzo ;t al.l I2010T) . In normal objects, the evi- 
dence of carbon is given by a weak absorption feature 
at ^6300 A attributed to C II A6580 (|Mazzalil [20011: 
Branch et aLll2003t iThomas et aLll2007H An alternative 
identification o f this a bsorption as Ha was suggested by 
iGaravini et"aT1 (|2004fl for SN 1999aa at 11 days before 
maximum light. The evidence, however, was not conclu- 
sive. Hydrogen from the companion star or a circum- 
stellar envelope could be mixed in the outer layers of 
the ejecta and possibly be detected in spectra at phases 
between 15 and 5 days befor e maximum. Based on spec- 
tral synthesis calculations, iLentz et al.1 (j2002[ ) showed 
that Ha can produce the 6300 A absorption at such 
early phases if they assumed the existence of ~10 -2 M Q 
of hydrogen at v > 15000 km s _1 . However, such a 
large amount of hydrogen is in contradiction with the 
limit of ~10~ 7 Mm imposed by mass-transfer models 
(jNomoto et al.ll2007l ). 

In this paper, we analyze the presence of the C II A6580 
line in a sample of pre- maximum sp ectra obtained by th e 
Carnegie Supernova Project (CSP: lHamuv et all [2006), 
and thereby infer the incidence of unburned material in 
normal SNe la. In Section [2] we make a brief description 
of the spectroscop i c data , most of which will be presented 
bv lFolatelli et al.l (|2011l ). and provide the evidence sup- 
porting the detection of C II lines. The alternative of hy- 
drogen is also studied. In Section [3] we analyze the spec- 
troscopic properties of SNe la with carbon in comparison 
with those which show no clear evidence of it. Section [4] 
presents a comparative study of photometric properties 
of SNe with and without carbon. Finally, in Section [5] we 
briefly discuss our results in light of the available explo- 
sion models. Note that betwee n the presentatio n of a pre- 
liminary report on this work (Folatclli 2010a ) and sub- 
mission of the current paper, iParrent et al.l (|201lD pub- 
lished similar results to the ones presented here, confirm- 
ing a significantly higher incidence of carbon in SNe la 
spectra than previously noted . In Sectional we com pare 
the findings of the paper by IParrent et al.l (|2011[ ) with 
the conclusions of our own study. 

2. DATA SET AND CARBON DETECTION 

We base the curr ent work on the sam ple of optical 
spectra presented bv lFolatelli et al.1 (|201lD . This sample 
consists of 588 spectra of 82 SNe la observed by the CSP 
between 2004 and 2009. In searching for signatures of 
unburned material, we only consider the pre-maximum 
phases. This is because any primordial C-0 rich mate- 
rial left by the explosion should lie in the outer parts 
of the ejecta, a nd is thus most eas ily detectable at the 
earliest stages (iTanaka et al.l 120081) . For comple t eness , 
we have added to the sample of rFolatelh eTaTl pOllh 
other spectra in the CSP database which were obtained 
before maximum light. These additional spectra belong 
to SNe 2008hj, 20091, and 2009cz. The present sample 
amounts to 157 spectra of 51 SNe la with ages between 
— 12 and days. Table [TJ provides the list of SNe used in 



this work. The same table also gives the epoch of the first 
spec trum , the class of carbon detection as defined in Sec- 
tion 12.11 and several photometric parameters which are 
used in Section [TJ Table [5] contains information of each 
spectrum. 

We have also included SN 2009dc which is a luminous, 
slowly expanding SN la suggested to arise from a 
super- C handrasekhar mass pr o genitor jYa manak a et all 
2009b : iTanaka et all \20m iSilverman et al.l 120111: 
Taubenberger et al.l 120111) . For consistency, we only 
include the pre-maximum spectra of this SN in Table [3J 
although it showed strong C II features which persisted 
after maximum light. In light of the peculiar nature of 
this object, we will treat it in Section [3] as a reference 
case to see whether there is a connection with the 
normal SNe la studied here. 

2.1. The Search for C II 

As mentioned in Section [I] carbon is the only element 
which can be identified with purely unburned material 
from a C+O WD progenitor. Oxygen instead could be 
due to a mixture of pri stine material and as hes from C 
burning. According to ITanaka et al.l (|2008l ) . the dom- 
inant carbon ionization state in the outer parts of the 
ejecta is C II. Therefore, in our study of pre-maximum 
spectra, we search for signatures of unburned material in 
the form of C II features. 

The four strongest C II lines at optical wavelengths are: 
C II A4267, A4745, A6580, and A7234. Of these, the most 
promine nt in SNe la spectra is expected to be C II A6580 
(e.g., see IHatano et al.lll999t rMazzalill200lD . Even if C II 
A4745 has a lower oscillator strength than C II A6580, 
it may become comparably strong at temperatures of 
«5000 K due to its lower excitation energy. However, 
at such low temperatures, the blue part of the spectrum 
is dominated by Ti II lines, which complicates the iden- 
tification of this carbon feature. For expansion velocities 
between 10000 and 20000 km s" 1 the C II A6580 line 
should produce absorption between 6170 A and 6370 A, 
which lies on the red side of the Si II A6355 line. We 
have therefore focused the search for unburned material 
in this wavelength region. 

Visual inspection of the spectra allowed us to divide 
the present SN sample into three main groups. First, 
those SNe which show no apparent absorption but rather 
a clear emission component of the Si II A6355 line. They 
are marked with "N" (i.e., "no detection") in Table [TJ 
Objects that do show a clear absorption feature, usu- 
ally centered at «6320 A are indicated with an "A" (i.e., 
"absorption" ) in Table [TJ An intermediate case is that of 
SNe which show no clear absorption but an apparently 
suppressed Si II A6355 emission component. Because of 
the flat emission profile, these are labeled "F" (i.e., "flat 
emission" ) in Table [TJ We note that the distinction be- 
tween groups A and F is not a quantitative one, and 
the specific classification may depend in each case on the 
quality of the available data. Nevertheless, this does not 
affect the results of the forthcoming analysis. Finally, 
for those cases which are inconclusive — typically due to 
data with low signal-to-noise ratio — we have added a 
"?" to the identification. In the cases of SNe 2007S, 
2008ar, 2008gp, and 2009ad, the available spectra do not 
cover the wavelength range under scrutiny and thus we 
are not able to classify them according to the presence of 
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the putative C II A6580 line. Figure Q] shows the spectra 
which belong to group A in the wavelength regions of 
the four C II lines. Figure [2] shows the spectra with flat 
Si II A6355 emission. In the left panel we include SNe in 
the F group while in the right panel objects are classified 
as "F?". As a comparison, several "non-detections", i.e 
SNe in group N, are shown in Figure [3J 

We have subsequently studied the spectral regions cor- 
responding to the other three C II lines in search of simi- 
lar features which may help consolidate the identification 
with this ion. Although in the case of some A-type ob- 
jects other absorptions at the right locations can be seen 
(Figure [IJ , these are usually very weak and it is difficult 
to clearly associate them with C II. In particular, the 
two lines on the blue part of the spectrum are subject to 
severe blending into stronger lines, mainly of Mg II, Si II, 
S II and Fe II. Moreover, an absorption matching the ex- 
pected position of the C II A4745 line is seen in most of 
the spectra, irrespective of the SN belonging to the A, F 
or N types. The absorption due to C II A7234, depending 
on the redshift and expansion velocity, may lie near the 
6900 A telluric band, which may affect the detection if 
the atmospheric feature is not accurately removed. 




4000 6000 7000 4000 6000 7000 

Rest frame wavelength [A] Rest frame wavelength [A] 



Fig. 1. — Spectra of all the SNe la classified in group A (with 
C II A6580 absorption) at the wavelength regions of the C II lines. 
The shaded bands mark the position of four C II lines blue-shifted 
between 9000 and 14000 km s" 1 . For clarity, the spectra are plot- 
ted as f v /v. Labels indicate the names of the SNe and their age in 
days relative to ZJ-band maximum light. The SNe are sorted from 
top to bottom in increasing order of Amis(B). The gaps in the 
spectra of SN 2009ab {right panel) are due to a separation between 
two detectors. 

A more careful identification of spectral lines can 
be carried out using the highly paramet erized super- 
nova synthetic spectrum code SYNOW (see lBranch et al.l 
l2002t and references therein). Examples of synthetic 
spectra compared with observations are shown in Fig- 
ures |4] through [T] Due to the approximations involved in 
SYNOW (spherical symmetry, sharp photosphere, LTE, 
Sobolev line formation, etc.), the synthetic spectra are 
not expected to produce a perfect match to the observed 
flux and line profiles. Nevertheless, SYNOW is useful 
for identifying spectral features with the species which 
produce them, specifically by reproducing the location 
of observed absorptions with predicted lines of a given 
ion. 

Table [3] provides a summary of the main parameters 
used to produce the synthetic spectra shown in Figures 0] 
throughUJ In order to simplify the analysis, several other 
parameters have been left fixed. In particular, we have 
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Fig. 2. — (Left) Spectra of all SNe la classified in group F (with 
flat Si II emission) at the wavelength regions of the C II lines. The 
shaded bands mark the position of four C II lines blue-shifted be- 
tween 9000 and 14000 km s — . For clarity, the spectra are plotted 
as f v /u. Labels indicate the names of the SNe and their age in 
days relative to j3-band maximum light. The SNe are sorted from 
top to bottom in increasing order of Ami5(B). (Right) The same 
for SNe la with F? classification. 
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Fig. 3. — (Left) Spectra of all SNe la classified in group N (with- 
out signs of C II A6580) at the wavelength regions of the C II lines. 
The shaded bands mark the position of four C II lines blue-shifted 
between 9000 and 14000 km s . For clarity, the spectra are plot- 
ted as fv/v. Labels indicate the names of the SNe and their age 
in days relative to S-band maximum light. The SNe are sorted in 
increasing order of Ami5(B). (Right) The same for SNe la with N 
classification and Si II expansion velocity greater than 13500 km 
s" 1 . In this case, the shaded bands mark the position of C II lines 
blue-shifted between 12000 and 20000 km s — 1 . The spectra are 
sorted from top to bottom in increasing order of Si II expansion 
velocity. 

adopted a power law of index 8 for the dependence of the 
optical depth on velocity, up to 30000 km s _1 . We have 
also fixed the excitation temperature to T oxc = 12000 K 
for all species. 

Figures S] and [S] show the cases of two objects in group 
A, namely SNe 2005el and 2009F at -7 and -5 days re- 
spectively. These SNe show different Ami5(B) values 
and correspondingly distinct spectroscopic properties. 
While SN 2005el is a normal event, SN 2009F presents 
fast-declining light curves an d belongs to the "Co ol" or 
"CL" type in the scheme of iBranch etaLl (|2006t) . The 
latter SN shows a relatively red spectrum and strong 
absorptions due to Ti II and O I, as compared with nor- 
mal SNe la. The synthetic spectra provide satisfactory 
matches to the observations, with small deviations which 
are not relevant to our analysis of carbon features. The 
ions producing the main spectral features are labeled (see 
Table|3]). In some cases, "detached" components of Ca II 
and Fe II have been included to reproduce high- velocity 
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(HV) absorption features (|Mazzali et al.ll2005l ). 

In the case of SN 2005el (see Figure HJ), C II can be 
invoked to reproduce two absorptions observed at about 
6300 and 7000 A. The most evident of the two is C II 
A6580 which appears sharper in the data than in the cal- 
culations. The assumed distribution of C II is slightly 
"detached" by 500 km s _1 above the photosphere in or- 
der to reproduce the location of the lines. The synthetic 
spectrum without C II (dotted line) fails to reproduce 
those two features. The identification of two absorptions 
provides extra support to the detection of carbon (see 
Section 12.41 for a comparison with other possible iden- 
tifications). In the blue, only C II A4745 produces a 
detectable absorption in the synthetic spectrum. This 
line roughly matches an observed absorption, although 
the identification is not clear due to the presence of S II 
and Fe III lines. Moreover, the synthetic spectrum in 
that region is dominated by an emission, which makes 
the identification even more doubtful. As will be shown 
below in the case of SN 2006ax, this emission could ten- 
tatively be suppressed by C III. No evident signature of 
C II A4267 is seen in either the observed or synthetic 
spectra at the expected position. 

We refer th e read er to a similar analysis done by 
iThomas et "all (|2007l ) on SN 2006D which is also in 
group A according to our classification. The authors dis- 
cuss whether non-LTE effects may modify the relative 
strength of C II lines to reproduce a strong and sharp 
absorption at about 4100 A which is not identifiable with 
any other species. 

The case of SN 2009F (see Figure [5} is similar to that 
of SN 2005el. As a consequence of SN 2009F being of 
the "Cool" type, the main differences in the synthetic 
calculations are a lower black-body temperature, and 
larger O I and Ti II optical depths. Additionally, no 
ion — including C II — is detached from the photosphere. 
Again the synthetic spectrum with C II provides a better 
match to the observations than the one without carbon, 
in particular by reproducing observed absorptions near 
6300 and 7000 A, although the apparent width of the 
C II A6580 line is too wide in the calculations. The iden- 
tification of C II lines in the blue part of the spectrum is 
even more difficult than in the case of SN 2005el due to 
the presence of strong Ti II lines. 

Figure [6] shows the case of an F spectrum of SN 2006ax 
obtained 10 days prior to maximum light. In this case, 
even if the absorption feature is not evident, a compari- 
son between the synthetic spectra with and without car- 
bon clearly shows that the presence of C II is a plausible 
explanation for the flat emission component of the Si II 
A6355 line. A weak absorption near 7000 A can be asso- 
ciated with C II A7234, as in the cases of SNe 2005el and 
2009F. We thus assume a common origin for the spectral 
properties of A and F spectra in the range around 6300 

A. 

Note also in Figure [6] that we have tentatively included 
C III in the synthetic spectrum. This adds an absorption 
at about 4500 A due to C III A4650 which serves to sup- 
press the emission components of Mg II and Si III lines 
to the blue. As mentioned above, C III can also be in- 
voked to improve the fit in this wavelength region in the 
case of SN 2005el, as shown in Figure 01 Although the 
presence of C III A4650 has been considered previously 



in the literature fsee | Garavini et al.ll2004t iBranch et al.l 
[20071 IParrent et al.ll2011lT its identification is only tenta- 
tive because it would require a very high temperature to 
excite this transition (29.5 eV). It should be noted that 
more elaborate models, for example those a pplied to the 
normal SN 2003du bv iTanaka et all (|2011f ). do not in- 
volve the presence of C III lines even when some carbon 
is included in the ejecta. 

For comparison, in Figure [7] we show the case of an N 
type spectrum. In this example, we show the spectrum 
of SN 2004eo obtained at —6 days along with a matching 
SYNOW spectrum. We see that, as opposed to the A and 
F cases, the match around 6300 A is good without the 
need to invoke C II. In a case like this we can confidently 
say that carbon is not detected. However, as shown in 
Section [2~2l in other cases classified as N or F?, carbon 
features may be hidden by noise or blending. 
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Fig. 4. — Spectrum of SN 2005el at —7 days (black solid line) 
which belongs to group A, and a matching SYNOW calculation 
(blue dashed line). The SYNOW synthetic spectrum without car- 
bon is shown for comparison near the C II A6580 and A7234 lines 
(green dotted line). The labels indicate the main ions which con- 
tribute to the formation of the observed features. The most impor- 
tant contributors for each feature are labeled closest to the spec- 
trum. Labels above the spectrum indicate that the ion was located 
"detached" from the photosphere. (Inset) Blow-up of the region 
around C II A6580 and A7234. The SYNOW spectrum containing 
H I instead of C II is also shown (red dotted line). The expected lo- 
cations of C II and Ha lines are indicated. While H I can reproduce 
the absorption at fs6300 A, C II produces a more solid identifica- 
tion by reproducing both absorptions at ss6300 and fs7000 A. 



2.2. Detection Statistics 

Assuming the identification of the 6300 A feature in A 
and F SNe is C II, we see a remarkably large incidence of 
unburned material in the current sample. If we consider 
the objects with spectroscopy obtained before or at —4 
days — the age of the oldest A spectrum — , from a total 
of 32 normal SNe la we have nine objects in group A 
and three in group F. This means 38% of the SNe are 
in the A or F groups, and 28% in A. Apart from the 
recent work by IParrent et al.l (|2011| ) who find a similar 
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Fig. 5.— Spectrum of the "Cool" SN 2009F at -5 days (black 
solid line) which belongs to group A, and a matching SYNOW 
calculation (blue dashed line). The SYNOW synthetic spectrum 
without carbon is shown for comparison near the C II A6580 and 
A7234 lines (green dotted line). See Figure [4] for a description of 
the labels. (Inset) Blow-up of the region around C II A6580 and 
A7234. The expected locations of C II lines are indicated. 
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Fig. 6. — Spectrum of SN 2006ax at —10 days (black solid line) 
which is classified in group F, and a matching SYNOW calculation 
(blue dashed line). The SYNOW synthetic spectrum without car- 
bon is shown for comparison near the C II A6580 and A7234 lines 
(green dotted line). See Figure [4] for a description of the labels. 
(Inset) Blow-up of the region around C II A6580 and A7234. The 
expected locations of C II lines are indicated. 

detection fraction of 30%, previous claim s for carbon in 
SNe la spectra have be en sporadic (see iThomas et al.l 
I2007t [Branch et aLll2007f ) . As shown below, we think this 
has to do with the relative scarcity of early-time spectra. 

Moreover, we note that the fractions derived above are 
lower limits because we count as non-detections several 
cases where the data are not suitable to fully discard 
the presence of carbon. Below we present three different 



Fig. 7. — Spectrum of SN 2004eo at —6 days (black solid line) 
which belongs to the N class, and a matching SYNOW calculation 
without carbon (blue dashed line). See Figure [4] for a description 
of the labels. (Inset) Blow-up of the region around C II A6580 and 
A7234. The expected locations of C II lines at about 1000 km s — 1 
above the Si II velocity are indicated. The same as above for the 
spectrum of SN 2004co at —6 days which belongs to the N class. 

effects which can work against the detection of C II lines: 
noise, expansion velocity, and age. 

2.2.1. Noise 

Clearly, noise in the data can hide weak features such 
as the ones under study here. In order to quantify the 
detectability of the absorption at 6300 A, we have per- 
formed equivalent width measurements. Strictly speak- 
ing, since the location of the actual continuum is un- 
known, these are pseudo-equivalent widths (pW). We 
simply trace a straight line along the absorption feature 
to mimic the continuum flux. An example of a pW mea- 
surement is shown in Figure [8] The pW is computed as 
the area of the absorption feature after normalizing the 
spectrum by a local pseudo-continuum. In cases when no 
absorption is seen, we derive an uncertainty in the pW 
value in the region around 6300 A in order to set a 3-ct 
upper limit. The measured values and upper limits are 
listed in Table [5J and shown in Figure GO 

We can see that, in general, pW > 1 A for objects 
in the A and F groups, and that most of the 3-cr upper 
limits of SNe with no C II detections lie below that value. 
However, in several cases the upper limits are comparable 
to or even larger than the measured pW of objects in 
groups A and F. This implies that some of the SNe la 
with carbon may have gone undetected because of noise 
in the spectra. 

2.2.2. Expansion Velocity 

The higher the expansion velocity, the wider the spread 
of spectral features in wavelength space, and thus the 
larger the impact of line blending. In particular, the hy- 
pothetical 6300 A absorption, which is much weaker than 
its neighboring Si II A6355 line, would become increas- 
ingly blended into the red absorption wing of the latter. 
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Fig. 8. — Example of a pscudo-cquivalcnt width measurement 
(pW) for the C II A6580 line. A pseudo-continuum is defined as a 
straight line traced between the flux peaks on each side of the ab- 
sorption feature (vertical lines). The spectrum is then normalized 
by this pseudo-continuum, and the pW is computed as the area of 
the normalized absorption (shaded region). 
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Rest — frame days relative to B maximum 

Fig. 9. — Pseudo-equivalent width of the absorption at «6300 
A which is attributed to C II A6580. Filled symbols are used for 
SNe la in A and F groups. Open triangles mark 3-cr upper limits 
for the rest of the spectra in the current sample. Open squares 
linked with lines show the results for synthetic spectra computed 
for different masses of car bon in the velocity range 10500 < v < 
15000 km s" 1 (see Section ROt. 

This would make the detection of carbon problematic 
when expansion velocities are high. 

To illustrate this problem, plotted in Figure [TU] are 
synthetic spectra obtained with SYNOW using different 
photospheric velocities. For comparison, data of differ- 
ent SNe with approximately matching expansion veloc- 
ities are shown along with the synthetic spectra. The 
top synthetic spectrum was tailored to match the ob- 
served one of SN 2006D at —6 days. We note that the 
C II A6580 absorption in the synthetic spectrum is wider 
than the observed one, but it reproduces well the ob- 



served intensity. This SYNOW spectrum is taken as the 
fiducial model to produce the rest of the synthetic spec- 
tra shown in the figure. In each case we have only varied 
the photospheric velocity and optical depth of Si II in or- 
der to match the location and strength of the Si II A6355 
absorption. Note that we have not attempted to produce 
an accurate match of each spectrum, but we have only 
tried to illustrate the effect of Doppler blending. 

It is clear from Figure [10] that, when the expansion ve- 
locity exceeds roughly 15000 km s _1 , a C II absorption 
as strong as or weaker than the one of SN 2006D would 
be very difficult to detect. In the case of SN 2004ef at —6 
days, the signal-to-noise ratio of the data makes the de- 
tection of the hypothesized C II line barely possible. For 
SN 20071e at —9 days, we can observe some "wiggles" on 
the Si II A6355 emission profile which may be compatible 
with a C II absorption. These weak features, along with 
the complex structure of the Si II absorption, disappear 
by the following observation of SN 20071e obtained at —4 
days. This led to our classification of this object as "F?" . 

The right-hand panel of Figure [3] shows several cases 
of non-detections of C II for SNe la with Si II velocities 
above 13500 km s _1 . In some of these cases, specifi- 
cally for SNe 2004dt and 2006X, the expansion velocities 
are so large that they would hide the weak C II absorp- 
tion observed in the case of type A SNe la. We note 
that the nature of SN-Ia explosions conspires with this 
effect against the detection of carbon. This is because ex- 
pansion velocities are the largest at the earliest phases, 
which implies that the line blending effect is worst at 
times when the plausible C II absorption is expected to 
be strongest (see Section 12.2.30 . 

As shown in Section [3j none of the spectra classified as 
A or F show a Si II A6355 velocity greater than «12500 
km s _1 . According to the brief analysis based on Fig- 
ure [TQl the detection of the C II A6580 line could be pos- 
sible for velocities above that value and up to «15000 km 
s _1 . This may indicate that the lack of carbon detections 
among spectra within this range of expansion velocities 
is not exclusively due to the effect of Doppler blending. 

2.2.3. Age 

All the claimed detections of C II in the spectra of nor- 
mal SNe la, as the o nes presented in this paper and in 
previous publications dThomas et aTll2007t iBranch et al.1 
l2007t [Parrent et aT1l2011[ ). correspond to early-time ob- 
servations, most commonly between one and two weeks 
before maximum light. This is mainly because carbon is 
expected to be present in the outer region of the ejecta 
which are visible at early phases. This is consistent with 
the temporal decrease in the strength of the absorption 
at «6300 A, as can be seen in Figure [9] The case of 
SN 2008bf, shown on the left panel of Figure [T] is a clear 
example of this evolution between —9 and —3 days. An 
additional spectrum of this SN, obtained at —1 day, is 
classified as N (see Tabled]). 

Sample statistics further support this picture. Fig- 
ure [TT] shows the fraction of carbon detections as a func- 
tion of age for the current sample. The fractions are 
simply the ratio of the number of SNe in group A (or 
A and F) to the total amount of objects whose spectra 
cover the region around 6300 A. As we go from —11 to 
— 1 day, the fraction decreases monotonically from about 
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Fig. 10. — The effect of Doppler blending on the possible de- 
tection of C II A6580. Four observed spectra are plotted (gray 
continuous lines) with the labels on the left side indicating the SN 
name and epoch in days relative to B maximum. SYNOW calcu- 
lations accompanying each spectrum show the cases with (dashed 
lines) and without (dotted lines) C II. The top synthetic spectrum 
was produced to fit the complete observed spectrum of SN 2006D 
at —6 days. This fiducial spectrum was used to produce the other 
synthetic spectra with different photospheric velocities. The labels 
on the right indicate this velocity in units of 10 3 km s — 1 . 

30% (40% for A+F) to zero detections. 

In the next section we show that some models predict 
a weaker absorption at ~6300 A before —11 days, which 
may lead to a lower detection fraction prior to that epoch. 
This highlights the importance of early observations in 
determining the incidence of unburned material in the 
ejecta of SNe la. 

2.3. The Mass of Carbon 

In an attempt to put constraints on the amount 
of unburned material which is necessary to explain 
the observed spectral features, we have computed syn- 
th etic spectra using the Monte Carlo co de developed 
bv iMazzali fc Lucvl (1 19931) and updated bv lLucvl (|19990 
and iMazzalil (|200(l ) . A recent description of the code 
and its appro ximations wa s give n in an analysis of 
SN 2003du bv ITanaka et aH (|2011f ). For that standard 
SN (jStanishev et al.ll2007| ) the authors find that a mass 
of carbon of M(C) = 6.8 x 1O~ 3 M (mass fraction 
X(C) = 0.002) in the velocity range 10500 < v < 15000 
km s" 1 is suitable to explain the suppressed emission 
peak of the Si II A6355 line observed at —11 days. More 
carbon is expected to be present in outer layers, but 
it would not produce an y noticeable e ffect in the spec- 
trum (see Figure CI of ' ITanaka etHI 120111) . The au- 
thors set an upper limit of 0.016 AT© for the total mass 
of carbon in the ejecta of SN 2003du. Based on this 
configuration we produce two additional models by only 
changing the mass of carbon in the velocity range given 
above to four times below and above the fiducial value of 
M{C) = 6.8 x 10 -3 M Q . The three models are employed 
to compute synthetic spectra at several phases between 
— 13 and +1 days from maximum light and thus evaluate 
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Fig. 11. — Fraction of SNe with C II as a function of age. Values 
show the ratio of SNe with spectra in the A (solid black) or A+F 
(shaded red) groups to the total number of SNe with spectra cover- 
ing the region around 6300 A. Horizontal error bars show the size 
of the epoch boxes used for computing statistics. Vertical error 
bars indicate the expected variation of the fraction for an assumed 
error of ±1 SN with carbon. 

the effect on the C II A6580 line. 

Figure [H shows the pW of C II A6580 as a function 
of SN phase for the synthetic spectra of different car- 
bon mass. As expected, the pW generally grows with 
mass. The range of observed pW is well represented by 
the assumed range of M(C). The lowest adopted mass 
(1.7 x 10~ 3 M©) would correspond to SNe with no car- 
bon detection, while the largest mass (2.7 x 1O~ 2 M0) 
roughly matches the strongest observed C II 6580 lines 
for SNe in group A at about one week before maximum 
light. 

In general, the models confirm the observed fact that 
the pW of carbon tends to decrease with phase. The 
exception to this is the increase in pW between the mod- 
eled spectra at —13 and —11 days. Close inspection of 
the synthetic spectra at —13 days shows a strong emis- 
sion component of the Si II line which washes away the 
weak absorption due to C II and produces a low pW as 
compared with that at —11 days. To a lower extent, the 
C II pW is additionally reduced by the fact that at —13 
days part of the carbon lies below the photosphere. Fur- 
thermore, a relatively larger expansion velocity at day 
— 13 enhances line blending, thus contributing to the de- 
crease in the measured pW . If this effect were generally 
valid, one would expect that the fraction of SNe with 
carbon detection should decrease at phases earlier than 
roughly —12 days. It should be noted, however, that 
the adopted model produces a worse match to the ob- 
served spectrum of SN 2 003du at —13 days t han at later 
phases (see Figure 2 of Tanak a et al.l |20TT1) . especially 
in the wavelength region under scrutiny here. Unfortu- 
nately our data do not cover such early phases and thus 
we cannot corroborate the results of the model. More 
early-time spectroscopy of SNe la is needed to settle this 
question. 

2.4. Hydrogen? 
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With the aid of SYNOW we have studied alternative 
identifications for the absorption around 6300 A in A- 
type SNe la. Other ions which may produce such absorp- 
tion at expansion velocities near the photospheric value 
are Ne I and H I. The strongest Ne I lines (at 6402 and 
6506 A) would only reproduce the absorption at «6300 
A if the expansion velocity were between 4000 and 9000 
km s — 1 . Such velocities are lower than those observed 
for SNe la before maximum light. 

On the other h and H I is a plausible alternative 
(|Lentz et al.l [2002). H-rich material may be present 
by mass transfer from a companion star in a single- 
degenerate progenitor system. If the observed absorption 
were due to H I this would rule out the double-degenerate 
scenario for a fraction of SNe la. 

As can be seen in the inset plot of Figure 2] (dotted 
lines), Ha at photospheric velocity is able to reproduce 
the 6300 A absorption as well as C II does. The weaker 
absorption at about 7000 A is, however, not reproduced. 
This makes the identification of H I less plausible than 
that of C II. With the given optical depth for Ha, the 
expected strength of H/3 is barely detectable, and it is 
not observed in the data. 

We have computed Monte Carlo synthetic spectra in- 
cluding H in order to model the observations. For this 
purpose we utilized the abundance distrib ution derived 
for SN 2003du at —11 days as given by iTanaka et al.l 
(poll . We added H with mass fraction up to 0.1 in the 
range of 10500 < v < 15000 km s _1 — i.e., right above 
the photosphere. We found that the absorption at 6300 
A for SNe in the A group can be reproduced with mass 
fractions X(H) ~ 0.05 which in this case translates into a 
total hydrogen mass of M(H) ~ 10 -2 M©. This is a very 
large amount of hydrogen which is incompatible with the 
theoretical picture of SNe la. In the single-degenerate 
scenario, models of progenitors with mass transfer pre- 
dict an upper limit of ~10~ 7 Mq whi ch is set by the 
ignition mass of H (jNomoto et al J 1200 71) . In the double- 
degenerate scenario, an even smaller amount of hydrogen 
is expected. 

Alternatively, a mass of the order of 0.01 M© of 
hydrogen could be stripped from the companion star 
in the single-degen erate scenario (|Marietta et al.l 120001: 
IPakmor et al.l [20081) . However, according to these simu- 
lations, the stripped hyd rogen is distributed in the inner 
layers (v < 4000 km s -1 ; IPakmor et a l. 2008), which are 
not visible in early-time spectra. 

Based on these observational facts and theoretical pre- 
dictions, even if the 6300 A absorption can be reproduced 
by hydrogen, we consider this possibility to be unlikely 
in favor of the C II identification. 

3. SPECTROSCOPIC PROPERTIES 

We now compare the spectral properties of the sam- 
ples of SNe with and without evidence of carbon. We 
first examine the expansion velocities. Figure [12] shows 
Si II A6355 expansion velocities as derived from the blue- 
shift of the absorption minimum. We see that the sample 
of normal SNe la with carbon appear to lie preferentially 
near the lower edge of the velocity distribution of the 
complete sample. Si II velocities are below wl2500 km 
s -1 in all of these case s. In terms of the classification 
by iBenetti et al.l ([2005) , none of the SNe with carbon 



are of high velocity-gradient (HVG) type. Accordingly, 
none of the SNe with car bon belong to the "Br oad Line" 
(BL) type as define d by iBranch et al.l (|2006l) (see also 
iFolatelli et "all 120 111) . As explained in Section HT2T21 the 
lack of carbon detection when velocities are high may 
be caused by an observational bias as a consequence of 
Doppler blending. However, the fact that in our sim- 
ulations C II lines are expected to be noticeable up to 
velocities of «15000 km s _1 , suggests that the prefer- 
ence for low velocities is real and not completely due to 
the observational bias. The possible implications of the 
velocity of carbon on SNe la models are briefly discussed 
in Section [5j 

Signatures of carbon are not exclusively found in nor- 
mal SNe la. In our sample, one fast-declining object, 
SN 2009F, belongs to group A (see Figure El). In the 
classification scheme of IBranch et al.l (|2006|) . this SN is 
of the "Cool" (CL) type. Three other CL SNe were ob- 
served at the range of epochs considered here, with no 
evidence of C II lines. The two SNe of the "Shallow Sili- 
con" (SS) type with early-time spectroscopy do not show 
evidence of carbon. 

In Figure Q~3] we present the expansion velocities of the 
bulk of carbon, measured from the absorption minimum 
located at «6300 A, under the assumption that it is due 
to C II A6580. On average, the evolution of C II A6580 
velocities is parallel to that of Si II A6355 velocities, and 
lying at «1000 km s _1 above. These velocities indicate 
that carbon appears slightly detached above the photo- 
sphere (see Section |2.1[) . If carbon were left unburned 
by a spherically symmetric explosion, then it would be 
expected to only b e present in the ou t ermost regions 
of th e ejecta (e.g., INomoto et al.l [19841 : llwamoto et al.l 
119991 predict a minimum velocity of 15000 km s _1 ). 
For current one-dimensional deflagration and delayed- 
detonation models, no carbon should be left at such low 
velocities. As discussed in Section [5] multi-dimensional 
effects may account for the presence of carbon deep in 
the ejecta. 

If the detection of carbon is related to prevalent asym- 
metries in the ejecta, we would expect to find evidence 
of other effects of asphericities. In addition to the dis- 
tinction between LVG and HVG SNe which was sug - 
gested as evidence of asymmetries ()Maeda et alj|2010"cf) . 
another possible effect would be t he presence of high - 
velocity components in the spectra (Tanaka et al. 2006). 
As was shown in Section 12.11 at the SN phases stud- 
ied here, high-velocity comp onents of Ca II and Fe II 
are commonly seen (see also IMazzali et aT] 12005). Fig- 
ure Q3] shows, however, that the presence or absence of 
high- velocity Ca II infrared (IR) triplet absorption is not 
obviously related to the presence of carbon. 

3.1. Super- Chandrasekhar SNe 

As a comparison we examine the CSP spectra of 
the peculiar SN 2009dc. This SN exhibits strong car- 
bon features in its pre-maximum spectra, which are 
also present until approximately a week after maxi- 
mum. From its spectral properties, slow light curves 
and exceedingly high luminosity, it has been suggested 
SN 2009dc resulted from a super-Chandrasekhar mass 
progenitor (lYamanaka et al.l l2009bl: ITanaka et al.l 120101 : 
iSilverman et alj|2011t iTaubenberger et al.ll2011l ). 
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Fig. 12. — Expansion velocities of Si II A6355 as a function of 
epoch relative to B maximum. Open black symbols mark SNe 
without car bon, with different sh ape for different spectral subtypes 
as given bv lFolatelli et al.l (f201l] V Filled symbols are used for SNe 
in the A and F groups. These symbols are the same as in Fig- 
ure [9] Data points for the same SN are linked with dotted lines. 
The shaded gray region marks the region within ltr about the aver- 
age ve locity evolution of normal SNe la as given bv lFolatelli et all 
(120111) . The yellow hashed region indicates the average and 1-a 
region for SNe with carbon. 
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Fig. 13. — Expansion velocities of C II A6580 {filled symbols) 
and Si II A6355 {open symbols) as a function of epoch relative to B 
maximum. The dashed line marks the average Si II A6355 velocity 
and the gray shaded area indicates the 1-tr dispersion. The solid 
line and hashed region mark the average and 1-cr dispersion of C II 
A6580 velocities. SN 2009dc is excluded from both averages. 

In Figure [T^] we compare the expansion velocities of 
Si II A6355 and C II A6580 derived from the spectra of 
SN 2009dc to our sample of normal SNe la. This figure 
illustrates SN 2009dc has lower Si II velocities than the 
normal objects by «2000 km s _1 . Interestingly, the rela- 
tion between C II and Si II velocities that was observed 
for normal SNe la also holds for SN 2009dc. In this 
case, the bulk of carbon also appears to be located above 
the silicon layer by «1000-2000 km s _1 . From a study 
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Fig. 14. — Spectra in the region of the Ca II IR triplet as a 
function of Doppler velocity assuming 8579 A as the effective rest 
wavelength. In the left panel, spectra classified in A and F groups. 
In the right panel, spectra of type N. The spectra are sorted accord- 
ing to their age (which is given in parentheses) . We see that cases 
with and without strong high-velocity components — near 20000 
km s — 1 — are present in both groups of SNe la. 

of a different sample of S Ne la, including a fe w super- 
Chandrasekhar candidates. iParrent et al.l (|2011l see their 
Figure 8) found a similar agreement of the relative loca- 
tion of carbon and silicon layers. This is an interesting 
finding because it may indicate that, even if the progeni- 
tors may be very different, the explosion produces similar 
relative distribution for the bulk of these elements. 

4. PHOTOMETRIC PROPERTIES 

This section presents a comparative analysis of the 
photometric properties of SNe la with and without sig- 
natures of carbon in the early-time spectra. For this 
purpose, a sample of 32 normal SNe la was selected 
with spectra obtained earlier than four days before max- 
imum light. Twelve of these objects belong to groups 
A and F, which are the ones considered to show car- 
bon in their ejecta. We b ase this study on the pho- 
tometric data publishe d by iContreras et al.l (|201Q[ ) and 
iStritzinger et al.l (1201 1|). The data wer e analysed with 
the SNooPy package ((Burns et al.l 120111 ) in order to de- 
rive light-curve parameters such as times and magnitudes 
at maximum light, and -B-band decline rates (Ami^(B); 
IPhillipslll99l . 

4.1. Decline rate and color 

Figure [15] (left panel) shows the distribution of light- 
curve decline rates, &mi§(B) for both groups. As a 
reference, we also show t h e dist ribution of the complete 
sample of iFolatelli et al.l (|2011l) . The group of SNe la 
with carbon covers a wide range of decline rates, includ- 
ing a fast-declining (1991bg-like) event. The average de- 
cline rate of all three samples (i.e., the complete sample, 
and those with and without carbon) is 1.2 ± 0.3 mag. 
Unfortunately, the samples sizes are too small to pur- 
sue a robust statistical test of the similarity among the 
distributions. 
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Shown in the right-hand panel of Figure [T5] is the dis- 
tribution of observed {B — V) pseudo-colors at maxi- 
mum light for SNe la with and without carbon, and 
for the reference sample. To avoid confusion with the 
intrinsically-red CL subclass of SNe la, we have excluded 
all SNe with Amis (2?) > 1-65 mag. We can see a ten- 
dency in the distribution of SNe with carbon to peak 
at a bluer color, although the statistics may be affected 
by the small number of cases. As poin t ed out recently 
(IPignata et all [2001 I Wang et al.l [2009t iFolev fc Kasenl 
120111: IMaeda et al.ll2011l) . SNe la with high expansion 
velocities tend to sh ow redder colors at maximum light. 
iTanaka et al.l (|2008l) suggested that the underlying cause 
of the difference in colors between HVG and LVG SNe 
is a variation in effective temperature which arises from 
a difference in photospheric velocity. The fact that SNe 
with carbon are all low-velocity objects may cause the 
observed difference in the color distributions. 

In order to further study possible differences in the 
intrinsic colors of SNe la with and without carbon, we 
have selected those object in each group which can be 
considered to have suffered little r eddening. Following 
the criteria of lFolatelli et al.l (|2010bf) , we picked SNe that 
(a) occurred in E/SO galaxies or that appeared outside 
the nuclei and arms of spiral galaxies, and (b) which did 
not show any sign of interstellar Na I D absorption in 
their early-time spectra. To these we added a few objects 
whose observed (B — V) colors one month past maximum 
light were compati ble with the in t rinsic c olor law given 
in Equation (2) of lFolatelli et"aH (|2010bD . The SNe in 
this low-reddening sample are indicated with a "Y" in 
column 8 of Table [Q 

Figure [16] shows the relation between color and decline 
rate for this sample of low-reddening SNe, with and with- 
out carbon absorption. Table 0] lists the resulting pa- 
rameters of the straight-line fits performed on different 
subsamples. We have adopted a Markov Chain Monte 
Carlo (MCMC) approach to perform the fits, with a full 
covariance matrix representation of the uncertainties in 
(B — V) and Amis(-B). We have also added an "intrinsic 
dispersion" term, 07, to the uncertainties which accounts 
for the extra dispersion in the data with respect to the 
straight-line fit. This term is left as a free parameter 
and its resulting values are given in Table [31 The "un- 
certainties" in 07 provide an estimate of the precision 
to which the MCMC method can determine this intrin- 
sic dispersion. Note that once again the Cool SNe have 
been excluded by requiring Amis(.B) < 1.6 mag. 

As shown in Table 0] and Figure [1^1 the relation be- 
tween color and decline rate for the combined sample of 
S Ne with low-reddening is compatible with the fit given 
by Folatcll i et al.l (|2010b| ). On the other hand, if the sam- 
ple is divided between those SNe with carbon and those 
without, a steeper relation is found for the former than 
for the latter SNe. Based on the posterior probability 
distributions resulting from the MCMC analysis, we find 
that the intercepts and slopes of both fits are different at 
the 80% and 90% confidence level, respectively. We must 
note that all of the SNe used in the fits — those which 
are assumed to be unreddened — happen to belong to the 
LVG class, so the difference cannot be explained by color 
variations between LVG and HVG SNe. Clearly, this 
difference must be confirmed using a larger SN sample. 
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Fig. 15. — {Left) Distribution of light-curve decline rates, 
Amis(B), for SNe la with carbon (black hashed) , without car- 
bon (l ight gray hashed), and the complete sample of lFolatelli et all 
(2011) (solid gray). (Right) Distribution of observed (B — V) 
pseudo-colors at maximum light for SNe with Ami^(B) < 1.6 mag. 
The shading code is the same as in the left panel. 




Am 15 (B) [mag] 

Fig. 16. — (B — V) pseudo-colors at maximum light versus 
Ami5(S). Solid symbols indicate SNe la with low reddening (see 
text). Circles mark SNe with carbon and squares, SNe without 
carbon. Straight-line fits to low-reddening objects are shown with 
a dashed line for SNe with carbon, dot-dashed line for SNe without 
carbon, and soli d line for both groups together. The dotted line is 
the fit given by Folatclli ct al. (2010b) using a different sample. 

4.2. Luminosity 

Now we compare the luminosities of SNe la with and 
without carbon. For this purpose, we compute distance 
moduli based on the redshifts and assuming a stan- 
da rd f?M~^A model of the universe (see Equation (5) 
of lFolatelli et al.l l2010bf ). which is accurate enough for 
the range of redshifts in the current sample. To avoid 
uncertainties larger than about 10% in the distances due 
to peculiar motions of galaxies in the local universe, the 
sample is cut by imposing that z > 0.01. 

Figure [T7] shows the histograms of B-band absolute 
peak magnitudes corrected for Milky- Way and host- 
galaxy extinction. In order to derive the latter cor- 
rection, color excesses E(B — V) were computed from 
the differences between observed (B — V) pseudo-colors 
at maximum light and t he law of intrinsic co lor ver- 
sus Amis (5) provided bv lFolatelli et al.l (|2010bD . Color 
excesses were converted into extinctions by multiply- 
ing by a total-to-selective absorption coefficient equiv- 
alent to the Galactic average of Rv =3.1. To avoid 
large uncertainties in the extinction correction, on the 
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left panel of Figure [T7] we cut the sample by requiring 
that (B — V) m ax < 0.2 mag, after correcting for Milky- 
Way extinction. We note an apparent difference in the 
extinction-corrected peak luminosities in the sense that 
SNe la with carbon are on average fainter than those 
without, and also than the reference sample. This dif- 
ference is however statistically marginal. The values are 
(Mb) = -19.25 ± 0.33 mag for 55 objects in the refer- 
ence sample, and (M%) = -19.23 ±0.38 mag for 15 SNe 
with no carbon detection. For the 6 SNe with carbon, 
this value is (M%) = -19.10 ± 0.21 mag, i.e., 0.13 mag 
fainter but well within the dispersion. 

A similar effect is seen in the histograms of the right- 
hand panel of Figure I17| where extinction-corrected, ab- 
solute peak magnitudes are shown, this time including 
objects with (B — V) max > 0.2 mag, but adopting a 
low value of the total-to-selective absorpti on coefficient of 
Rr = 2 .74, as derived from Hubble fits bv lFolatelli et al.l 
(|2010bl ). The differences in luminosity may in part be 
due to inaccurate extinction corrections caused by a dif- 
ference in intrinsic colors as that seen in Figure 1161 If, 
as suggested by Figure [THl SNe with carbon show bluer 
intrinsic colors then the extinction corrections for this 
group would be underestimated. In fact, if we adopt 
a different extinction correction for each group of SNe 
based on the intrinsic-color laws given in Table HJ the 
differences in average extinction-corrected peak magni- 
tudes nearly vanish. This may indicate that the differ- 
ences between both groups are related to differences in 
intrinsic colors rather than intrinsic luminosities. 

In Figure [15] we show the distributions of Hubble resid- 
uals in B band, AMb- For this p urpose, we ado p ted the 
parameters of fit 2 in Table 8 of iFolatelli et al.l ([2010bD 
which correspond to a fit of distance moduli versus 
Amis(B) and (B — V) pseudo-colors. The discrepancy 
between SNe with and without carbon in this case is seen 
in the same direction as above, although reduced to 0.08 
mag. The mean residuals and their standard deviations 
are (AMb) = 0.03±0.18mag and (AM B ) = -0.05±0.21 
mag for SNe la with and without carbon, respectively. 

Finally, we consider the effect of reddening on absolute 
peak magnitudes. For this purpose, we take the values 
of My at maximum light corrected for decline rate by 
adoptin g the slope of fit 7 in Table 9 of IFolatelli et al.l 
(2010b). In Figure [H we show these corrected peak 
magnitudes versus E(B — V) color excesses derived as 
explained above. With the current sample, we find no ev- 
idence of different extinction properties between SNe la 
with and without carbon. The data favors a low value of 
Rv ~ 2.1 as has been found in previous studies of SNe la 
data. 

5. DISCUSSION & CONCLUSIONS 

The most striking result of this paper is the large inci- 
dence of carbon lines in SNe la. In the CSP sample, we 
have found evidence of C II A6580 absorption in at least 
30 % of the objec t s, whi ch agrees closely with the findings 
of iParrent et alJ (|2011l ) based on a completely indepen- 
dent sample of publicly-available spectra. The identifi- 
cation of this absorption with other elements, such as 
hydrogen, is considered unlikely (see Section 12^4")) . In ad- 
dition, we have shown that nearly 10% of the SNe la in 
the CSP sample show a suppressed emission component 
of the Si II A6355 line, which SYNOW models indicate 




Fig. 17. — Distributions of extinction-corrected B-band abso- 
lute peak magnitudes. {Left) Extinction corrections derived using 
Ry = 3.1 and cutting the sample so that (B — V) < 0.2 mag. 
(Right) Including the complet e range of ( B — V) col ors, but adopt- 
ing a low value of Rp = 2.74 (Folatelli et al. 2010bh. The shading 
code is the same as in Figure fl5l 
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Fig. 18. — Di stribution of B-band H ubble residuals based on Fit 
2 of Table 8 in IFolatelli et al.l 12010bfl . Shading code is the same 
as in Figure [TS] 

may be caused by absorption due to the same C II A6580 
line. Taken together, these results imply that nearly half 
of all SNe la show carbon in their spectra at early phases, 
a finding that contradicts previous con clusions (e.g., see 
iThomas et all [20071 [Branch et al. 1 120071) . The CSP spec- 
tra show clearly that the detection of C II absorption 
is a strong function of the age, falling from >40% at 10 
days before B maximum to zero by the time of maximum 
itself. These results almost certainly represent a lower 
limit on the presence of unburned material in SNe la since 
very few events have been observed at even earlier phases. 
Moreover, low signal-to-noise observations, or line blend- 
ing when expansion velocities are high, also contribute 
to the difficulty of detecting the C II A6580 line. As we 
have shown, the latter effect can explain at least in part 
the absence of detectable car bon absorption in m embers 
of the high velocity gradient (|Benetti et al.ll2005D or BL 
([Branch et al J 1200a ) subclasses found bv IParrent et all 
((201 If ) and ourselves. As additional data are obtained at 
epochs of —10 days or earlier, it will be possible to cor- 
roborate whether, as we suspect, the majority of SNe la 
show carbon at some stage of their evolution. The theo- 
retical implications of such findings are potentially large, 
in particular as pertains to the mechanism by which the 
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Fig. 19. — Decline-rate corrected V-band absolute peak magni- 
tude versus E(B — V) color excess for SNe la with (circles) and 
without (squares) carbon. The solid line is a straight-line fit to 
all the data points, while the dashed and dot-dashed lines corre- 
spond to the SNe la with and without carbon, respectively. As a 
reference, the dotted line corresponds to the case of Ry =3.1. All 
three fits are in agreement with each other and favor a low value 
of Ry. Three SNe with the largest reddening — those which weigh 
most in the fits — are indicated. 

explosive name propagates inside the WD star. 

The second remarkable fi nding from ou r stud y of the 
CSP spectra and the work of lParrent et al.l (|2011[) is that, 
assuming that the absorption observed at about 6300 A 
is due to C II A6580, then carbon must be present in very 
deep regions, corresponding to velocities as low as v w 
11000 km s" 1 . whic h is right above the bulk of silicon 
(jMazzali et al.l 120071 ). This is well below the expected 
limit imposed by one-dimensional models, and points di- 
rectly to large mixing effects and/or possible departures 
from spherical symmetry or dumpiness. 

In spherically symmetric models, irrespective of the 
details of the flame propagation (deflagration or detona- 
tion), the production of a lar ge amount of 56 Ni as t ypi- 
cally observed (~ 0.6M e ; see iStritzinger et aLll2006l ) re- 
quires that the strength of the flame — either its speed or 
the detonation transition density — is high, which must 
lead to the total consumption of carbon below 1500 
km s- 1 (e.g. INomoto et al.|[l98l llwamoto et al.l 119991 ). 
and even below 20000 k m s _1 in detonation models (e.g. 
IShigevama et al.l [199 2). The situation is different for 
non-spherical models. For examp le, the of f-center igni- 
tion m odel of lMaeda et al.l ([2010bf ) (see also lKasen et al.l 
120091) synthesizes 0.54 M of 56 Ni, despite the relatively 
low density at which the detonation is triggered. The low 
transition density in their model results in less burning in 
the outer layers, and thereby the existence of carbon with 
a mass fraction of ~ 0.1 at velocities as low as ^13000 
km s . This is still larger than the observed velocity, 
but suggests that the non-spherical effects may be im- 
portant in understanding the detection of carbon deep 
in the ejecta. This argument is based on the effect of 
a global asymmetry produced by the off-center ignition. 
Similar conclusions may be reached from possible effects 
of mixing or clump formation in multi-dimensional mod- 



els. 

The observed rate of carbon detections could thus pro- 
vide information on the distribution of C-rich material, 
and the amount and size of the clumps. Unfortunately, 
the current uncertainties in the detection rate complicate 
these geometrical considerations. Nevertheless, if car- 
bon were present in a single clump or asymmetric blob, 
the tight dispersion in expansion velocities of C II A6580 
compared with sili con observed in Fig ure 1131 would be 
surprising (see also lParrent et al.ll201lD . 

In view of recent results which connect asymmetries 
in the ejecta with the observe d spec t roscopic and photo - 
metric diversity (jMaund et al.ll2010t iMaeda et al.ll2011l) . 
it is worth studying possible connections with the pres- 
ence of unburned material in the ejecta. Breaking our 
sample of CSP SNe la into two groups depending on 
the presence or absence of high- velocity components of 
Ca II absorption — a possible indicator of asymmetry 
(|Tanaka et al.l [2006]) — we find no significant distinction 
in the frequency of detectable carbon. There is an indi- 
cation in our data that SNe la with carbon show slightly 
bluer colors and lower luminosities at maximum light 
when compared with the rest of the sample. These are 
interesting hints of a possible origin of the observed di- 
versity arising from incomplete burn ing which may be re - 
lated to the asymmetry scenario of lMaeda et all poTfl) . 
but the statistical significance of these results is limited 
by the sample size. 

As shown in Section 12.31 modeling indicates that low 
masses of carbon in the range M(C) ~ 10 -3 - 10 -2 M© 
are enough to reproduce the observed spectral features, 
although high-velocity components, undetectable due to 
blending with the Si II A6355 line, could increase this 
value by as much as an order of magnitude (Tana ka et al.l 
120111) . Such small amounts of unburned material would 
make it difficult to detect other possible signatures in the 
observations. For instance, the non-detecti on of C I fea- 
tures in near-infrared spectra reported by iMarion et al.l 
(2006) could be expl ained by abundance and ionization 
state considerations (jTanaka et al. 2008). 

Another interesting possibility is an enhancement of 
oxygen versus magnesium in SNe which suffered less com- 
plete burning. While oxygen may in part be primordial, 
magnesium is only a product of carbon burning. De- 
termining the distribution of these elements in the outer 
eject a may shed light on t he evolution of the burning pro- 
cess (jMazzali et al.l [20081) . This could be attained with 
a temporal followup of the dominant spectral lines of 
these species, which at early times are those of O I and 
Mg II. Unfortunately, the study of such lines is compli- 
cated from an observational viewpoint. The main O I line 
at optical wavelengths usually coincides with the telluric 
A band, which makes it difficult to obtain an accurate 
measure of its strength from the ground. Other O I and 
Mg II lines appear blended with each other or with lines 
of other elements (e.g, see Figure 01 . I R this sense, the 
study of near-IR spectra may offer the poss ibility of find- 
ing is olated oxygen and magnesium lines ([Marion et al.l 
[20091 ). 

Figure Qj] clearly shows that the definitive assessment 
of the incidence of unburned material in SNe la require 
a survey of early-time spectra — obtained at least ten 
days before maximum light — covering the region of the 
C II lines. High signal-to-noise data are required in or- 
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der to explore the possible presence ol carbon in high- 
velocity SNe la, in particular by observing the relatively 
unblended but weak C II A7234 line. 
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—10.5 


F 


15 


002(028) 


-0.048(028) 


1 


02 


02) 


—0.022(057) 


Y 


0.01674 


-19 


40 


13) 


—0.11(15) 


2006dd 


-10.8 


A 


12 


241(009) 


-0.071(013) 


1 


08 


01) 


-0.053(052) 


Y 


0.00587 


-19 


01 


05) 


0.29(06) 


2006gt 


-0.3 


N 


18 


229(020) 


0.256(024) 


1 


85 


29) 






0.04477 


-18 


15 


06) 


-0.31(23) 


2006hx 


-9.6 


N 


17 


510(031) 


-0.157(073) 


1 


15 


04) 


-0.147(088) 




0.04549 


-18 


90 


06) 


0.58(21) 


2006kf 


-5.1 


N 


15 


817(005) 


-0.016(033) 


1 


58 


01) 


-0.058(060) 


Y 


0.02130 


-18 


91 


11) 


-0.13(14) 


2006mr 


-2.3 


N 


15 


345(011) 


0.708(121) 


1 


78 


04) 






0.00587 


-15 


90 


05) 


0.75(34) 


2006os 


— 1.9 


N 


17 


500(020) 


0.314(024) 


1 


14 


07) 


0.289(055) 




0.03281 


-18 


19 


07) 


-0.21(11) 


2007A 


-6.2 


F 


15 


695(012) 


0.222(033) 





92 


02) 


0.259(060) 




0.01765 


-18 


52 


13) 


0.09(16) 


2007S 


— 11.9 




15 


789(013) 


0.371(031) 





95 


03) 


0.405(059) 




0.01388 


-18 


22 


15) 


-0.03(17) 


2007af 


— 11.1 


A 


13 


330(014) 


0.196(030) 


1 


20 


03) 


0.200(059) 




0.00546 










2007ai 


-2.7 


N 


16 


805(009) 


0.119(011) 





74 


01) 


0.178(051) 




0.03166 


-18 


87 


07) 


0.16(08) 


2007as 


-2.2 


N 


15 


387(006) 


0.052(008) 


1 


22 


01) 


0.054(051) 




0.01757 


-19 


01 


12) 


-0.14(12) 


2007ax 


-2.1 


N 


16 


201(029) 


0.684(031) 


1 


96 


09) 






0.00686 










2007bd 


-7.5 


N 


16 


551(022) 


0.029(024) 


1 


23 


06) 


0.029(056) 




0.03102 


-19 


12 


07) 


-0.21(11) 


2007jg 


-0.2 


N 


17 


181(045) 


0.028(048) 


1 


31 


08) 


0.019(070) 




0.03713 


-18 


80 


08) 


0.06(16) 


20071c 


-9.2 


F? 


13 


875(009) 


0.292(011) 





97 


07) 


0.323(051) 




0.00672 










2007on 


-3.9 


F? 


13 


023(008) 


0.088(011) 


1 


86 


03) 




Y 


0.00649 










2008ar 


-8.1 




16 


259(038) 


0.100(054) 





88 


07) 


0.142(073) 




0.02615 


-19 


06 


09) 


-0.09(18) 


2008bc 


-10.1 


N 


14 


624(007) 


-0.038(015) 





82 


02) 


0.011(052) 




0.01509 


-19 


18 


14) 


-0.08(15) 


2008bf 


-8.5 


A 


15 


779(022) 


-0.066(031) 





87 


04) 


-0.023(059) 


Y 


0.02403 


-19 


36 


09) 


0.08(13) 


2008bq 


0.3 


N 


16 


661(006) 


0.052(008) 





92 


01) 


0.090(051) 




0.03400 


-19 


18 


06) 


-0.10(07) 


2008fp 


-3.7 


A 


13 


829(020) 


0.491(028) 





84 


05) 


0.538(057) 




0.00566 










2008gl 


-1.5 


N 


16 


845(022) 


0.050(031) 


1 


29 


04) 


0.043(059) 




0.03402 


-18 


90 


07) 


-0.09(12) 


2008gp 


-6.0 




16 


455(031) 


0.033(043) 





97 


11) 


0.065(066) 




0.03341 


-19 


28 


08) 


-0.19(16) 


zUUotlJ 


— 0.0 


iN 




o44^Uz / ) 


U.UUo^Uoo ) 


1 


04 


05) 


U.UZo^UOo J 




U.Uo/oy 


-19 


14 


07) 


— U.Uo^io ) 


2008hv 


-6.6 


A 


11 


799(006) 


0.067(009) 


1 


25 


01) 


0.065(051) 


Y 


0.01255 


-18 


99 


16) 


-0.19(16) 


2009D 


-3.9 


N? 


15 


789(016) 


0.018(023) 





84 


03) 


0.065(055) 


Y 


0.02501 


-19 


31 


09) 


-0.08(11) 


2009F 


-5.4 


A 


16 


929(025) 


0.635(036) 


1 


97 


04) 




Y 


0.01296 


-16 


75 


17) 


-0.04(20) 


20091 


-5.1 


F 


18 


247(081) 


0.702(115) 


1 


14 


16) 


0.713(125) 




0.02617 


-16 


92 


12) 


0.22(36) 


2009Y 


-5.5 


N? 


14 


042(032) 


0.114(045) 





88 


06) 


0.156(067) 




0.00935 


-19 


09 


22) 


-0.16(26) 


2009aa 


-7.8 


F? 


16 


355(023) 


0.030(033) 


1 


12 


04) 


0.044(060) 


Y 


0.02731 


-19 


06 


08) 


-0.07(12) 


2009ab 


-10.2 


A 


11 


652(046) 


0.065(065) 


1 


21 


09) 


0.068(082) 




0.01117 


-18 


65 


21) 


0.19(28) 


2009ad 


-8.3 




16 


074(096) 


-0.006(136) 





88 


17) 


0.036(145) 




0.02840 


-19 


33 


12) 


-0.07(41) 


2009ag 


-3.5 


F? 


11 


589(041) 


0.241(058) 





97 


08) 


0.273(077) 




0.00864 










2009cz 


-4.1 


N 


15 


786(053) 


0.073(075) 


1 


05 


10) 


0.095(090) 




0.02114 


-19 


06 


11) 


-0.13(25) 


2009dc 


-7.7 


SC 


15 


141(066) 


0.145(093) 





59 


11) 


0.222(106) 




0.02139 


-19 


68 


12) 


-0.61(29) 



Note. — Columns: (1) SN name; (2) Epoch of the first spectroscopic observation relative to S-band maximum light; (3) C 11 detection 
type (sec Section |2. Absorption (A); Flat emission (F); No detection (N); (4) MW-extinction corrected, corrected, apparent B- 
band peak magnitude (uncertainties in thousandth of mag); (5) Observed pseudocolor at maximum light, corrected for MW reddening 
(uncertainties in thousandth of mag); (6) Observed decline rate, Amis(B) (uncertainties in hundredth of mag); (7) Host galaxy color 
excess (uncertainties in thousandth of mag); (8) SNe with "Y" belong to the low-reddening sample; (9) Heliocentric rcdshift from the 
NASA/IPAC Extragalactic Database (NED); (10) S-band absolute peak magnitude (uncertainties in hundredth of mag); (11) B-band 
Hubble residual (sec Scction[4] uncertainties in hundredth of mag). 
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TABLE 2 
List of pre-maximum spectra 



UT Date 


MJD 


Phase 


Wavelength 


Fiesol. 


C II 


pPF(6300) 


Rcf. 








Ran trr> I A 1 




Class 


(A) 




l'\ "> 
( L ) 








C5~l 


y°) 


(71 

(') 










SN 2004dt 










2004-08-12 


53230.22 


-8.7 


4320 - 10315 




N 


< 1.11 


2 


2004-08-12 


53230.67 


-8.2 


3140 - 8589 




N 


< 0.75 


2 


2004-08-13 


53231.23 


-7.7 


3946 - 9702 




N 


< 1.02 


2 


2004-08-13 


53231.63 


-7.3 


3170 - 10250 




N 


< 0.66 


2 


2004-08-15 


53233.64 


-5.3 


3142 - 10400 




N 


< 0.33 


2 


2004-08-15 


53233.69 


-5.3 


3617 - 8879 




N 


< 0.39 


2 


2004-08-16 


53233.71 


-5.2 


3087 - 10211 




N 


< 0.24 


2 


2004-08-16 


53234.57 


-4.4 


3236 - 10678 




N 


< 0.45 


2 


2004-08-18 


53235.69 


-3.3 


3334 - 8878 




N 




2 


2004-08-18 


53236.54 


-2.5 


3217 - 10144 




N 




2 


2004-08-19 


53237.18 


-1.8 


3317 - 8894 




N 




2 


2004-08-20 


53238.22 


-0.8 


3318 - 8822 




N 




2 


2004-08-21 


53239.23 


0.2 


3337 - 8822 




N 




2 








SN 2004ef 










2004-09-07f 


53255.66 


-8.5 


3780 - 7323 


5.0 


N 


< 2.01 


1 


2004-09-08 


53256.68 


-7.5 


4000 - 7183 


5.0 


N 


< 0.99 


1 


2004-09-09f 


53257.77 


-6.5 


3650 - 6850 


14.0 


N 


< 1.05 


1 


2004-09-10 


53258.68 


-5.6 


4000 - 7114 


5.0 


N 


< 8.73 


1 


2004-09-14 


53262.61 


-1.8 


3800 - 9235 


8.0 


N 




1 








SN 2004eo 










2004-09-19 


53267.52 


-11.1 


4530 - 9618 


8.0 


N 


< 1.47 


1 


2004-09-19 


53268.05 


-10.6 


3501 - 9204 


0.8 


N 


< 0.30 


3 


2004-09-20 


53268.51 


-10.1 


3800 - 9235 


8.0 


N 


< 0.42 


1 


2004-09-24 


53272.74 


-5.9 


3341 - 10700 


0.2 


N 


< 0.21 


3 


2004-09-27 


53276.39 


-2.3 


3067 - 10676 


0.6 


N 




3 








SN 2004ey 










2004-10-25 


53303.55 


-0.8 


3600 - 9000 


14.0 


N 




1 








SN 2004gs 










2004-12-13 


53352.76 


-3.3 


3800 - 9235 


8.0 


N 




1 








SN 2004gu 










2004-12-19 


53358.83 


-2.9 


3600 - 9000 


14.0 


N 




1 








SN 2005M 










2005-02-04 


53405.71 


-0.1 


3800 - 9235 


8.0 


N 




1 








SN 2005W 










2005-02-03 


53404.51 


-7.5 


5760 - 7380 


5.0 


F? 


< 2.52 


7 








SN 2005ag 










2005-02-12 


53413.89 


0.3 


3800 - 9235 


8.0 


N? 




1 








SN 2005bl 










2005-04-16 


53475.64 


-8.4 


3800 - 9235 


8.0 


F? 


< 2.70 


1 


2005-04-17 


53477.60 


-6.5 


3450 - 8803 


10.0 


F? 


< 4.98 


4 


2005-04-19 


53479.63 


-4.5 


3800 - 9235 


8.0 


F? 


< 5.97 


1 


2005-04-19 


53479.70 


-4.4 


4960 - 10449 


11.0 


F? 


< 0.87 


4 








SN 2005bo 










2005-04-18 


53478.73 


-1.0 


3800 - 9235 


8.0 


N 




1 


2005-04-19 


53479.73 


-0.0 


3800 - 9235 


8.0 


F? 




1 








SN 2005el 










2005-09-26 


53639.82 


-7.3 


3800 - 9235 


8.0 


A 


3.58 ±0.26 


1 


2005-09-27 


53640.84 


-6.3 


3800 - 9235 


8.0 


A 


3.44 ±0.18 


1 








SN 2005eq 










2005-10-05 


53648.81 


-5.2 


3800 - 9235 


8.0 


N 


< 1.59 


1 








SN 2005hc 










2005-10-18 


53661.71 


-5.2 


4000 - 10200 


9.0 


N? 


< 7.83 


1 


2005-10-18 


53661.71 


-5.2 


5800 - 10200 


9.0 


_ 




1 








SN 2005ke 










2005-11-23 


53697.72 


-1.4 


3780 - 7289 


5.0 


N 




1 


2005-ll-24f 


53698.72 


-0.4 


3780 - 7290 


5.0 


N 




1 


2005-11-25 


53699.61 


0.4 


3200 - 9589 


14.0 


N 




1 








SN 2005ki 










2005-11-23 


53697.85 


-7.6 


3780 - 7289 


5.0 


F? 


< 2.40 


1 


2005-ll-25f 


53699.82 


-5.7 


3780 - 7289 


5.0 


N 


< 0.72 


1 






SN 2006D 










2006-01-16 


53751.83 


-6.1 


3200 - 5300 


6.0 


- 




1 


2006-01-16 


53751.84 


-6.1 


4000 - 10200 


9.0 


A 


3.74 ±0.21 


1 


2006-01-16 


53751.85 


-6.0 


5800 - 10200 


9.0 


A 


3.34 ±0.22 


1 








SN 2006X 










2006-02-09 


53775.71 


-10.4 


4135 - 6795 




N 




5 


2006-02-09 


53775.81 


-10.3 


3834 - 8139 




N 




5 


2006-02-11 


53777.65 


-8.5 


4175 - 6794 




N 




5 


2006-02-12 


53778.71 


-7.4 


4171 - 6794 




N 




5 


2006-02-13 


53779.63 


-6.5 


4178 - 6795 




N 




5 


2006-02-13 


53779.84 


-6.3 


3300 - 5300 


6.0 


N 




1 


2006-02-13 


53779.86 


-6.3 


4000 - 10200 


9.0 


N 


< 0.48 


1 
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TABLE 2 — Continued 



UT Date 
V 1 ! 


MJD 


Phase 
\ a ) 


Wavelength 

Ran era { A ^ 
-TV,ailgC \r\) 


Resol. 

1 k\ 
( A ) 

a\ 
\°) 


C II 

Place 


pW(6300) 
C71 


Rcf. 


2006-02-13 


53779.86 


-6.3 


5800 - 10200 


9.0 


N 


< 0.84 


1 


2006-02-18 


53784.59 


-1.5 


3902 - 8182 




N 




5 


2006-02-19 


53785.77 


-0.4 


4182 - 6800 




N 




5 








SN 2006ax 










2006-03-22f 


53816.76 


-10.5 


3800 - 9235 


8.0 


F 


1.21 ±0.47 


1 


2006-03-23 


53817.65 


-9.6 


3800 - 9235 


8.0 


F 


1.29 ±0.26 


1 


2006-03-24 


53818.75 


-8.5 


3800 - 9235 


8.0 


F 


1.91 ±0.35 


1 


2006-03-30 


53824.78 


-2.6 


3800 - 9235 


8.0 


N 


< 0.54 


1 


2006-03-31 


53825.73 


-1.6 


3800 - 9235 


8.0 


F? 




1 


2006-04-02 


53827.64 


0.2 


3800 - 9235 


8.0 


N 




1 








SN 2006dd 










2006-06-21 


53907.92 


-10.8 


3413 - 9610 


3.0 


N 




6 








SN 2006gt 










2006-09-25 


54003.75 


-0.3 


3800 - 9235 


8.0 


F? 




1 








SN 2006hx 










2006-10-05 


54013.65 


-9.6 


4000 - 10200 


9.0 


N 


< 1.05 


1 


2006-10-10 


54018.77 


-4.7 


3800 - 9235 


8.0 


N 


< 0.99 


1 








SN 2006kf 










2006-10-27 


54035.81 


-5.1 


3800 - 9235 


8.0 


N 


< 1.26 


1 








SN 2006mr 










2006-11-09 


54048.65 


-2.3 


3785 - 6131 


2.0 


N 




6 


2006-11-09 


54048.66 


-2.3 


5675 - 9976 


4.0 


N 




6 








SN 2006os 










2006-11-22 


54061.68 


-1.9 


3800 - 9235 


8.0 


N 




1 








SN 2007A 










2007-01-06 


54106.54 


-6.2 


3500 - 5300 


6.0 


— 




1 


2007-01-06 


54106.56 


-6.2 


4000 - 10200 


9.0 


F 


0.29 ±0.53 


1 


2007-01-06 


54106.57 


-6.2 


5800 - 10200 


9.0 


F 


1.85 ±0.64 


1 








SN 2007S 










2007-01-31 


54131.82 


-11.9 


3500 - 5300 


6.0 


_ 




1 


2007-02-12 


54143.73 


-0.1 


3800 - 9235 


8.0 


N 




1 








SN 2007af 










2007-03-04 


54163.82 


-11.1 


3300 - 5300 


6.0 


_ 




1 


2007-03-04 


54163.84 


-11.1 


4000 - 10200 


9.0 


A 


1.33 ±0.09 


1 


2007-03-04 


54163.86 


-11.0 


5800 - 10200 


9.0 


A 


0.99 ±0.11 


1 


2007-03-10 


54169.78 


-5.2 


3789 - 10889 


6.0 


A 


0.50 ±0.03 


1 


2007-03-14 


54173.80 


-1.2 


3429 - 9645 


8.0 


N 




1 


2007-03-14 


54173.81 


-1.1 


4962 - 9645 


8.0 


N 




1 








SN 2007ai 










2007-03-10 


54169.88 


-2.7 


4528 - 10187 


6.0 


N 




1 


2007-03-12 


54171.89 


-0.7 


4528 - 10192 


6.0 


N 




1 








SN 2007as 










2007-03-19 


54178.57 


-2.2 


3419 - 9626 


8.0 


N 




1 








SN 2007ax 










2007-03-26f 


54185.56 


-2.1 


3400 - 9608 


8.0 


F? 




1 






SN 2007bd 










2007-04-09 


54199.64 


-7.5 


3502 - 9697 


4.0 


N 




1 


2007-04-12 


54202.55 


-4.6 


3800 - 9235 


8.0 


N 


< 0.99 


1 


2007-04-17 


54207.59 


0.2 


3800 - 9235 


8.0 


N 




1 








SN 2007bc 










2007-04-09 


54199.62 


-0.5 


3506 - 9694 


4.0 


N 




1 








SN 2007jg 










2007-09-22 


54365.88 


-0.2 


3719 - 9961 


5.0 


N 




1 








SN 20071e 










2007-10-16 


54389.65 


-9.2 


3412 - 9623 


8.0 


F? 


< 0.39 


1 


2007-10-21 


54394.58 


-4.3 


3352 - 9560 


8.0 


N 


< 0.60 


1 


2007-10-21 


54394.60 


-4.3 


4704 - 9561 


8.0 


N 


< 0.54 


1 








SN 2007on 










2007-11-11 


54415.79 


-3.9 


3800 - 9235 


8.0 


F? 


< 0.24 


1 


2007-11-14 


54418.84 


-0.8 


3800 - 9235 


8.0 


N 




1 








SN 2008ar 










2008-02-29 


54525.87 


-8.1 


3743 - 10730 


6.0 


N 




1 








SN 2008bc 










2008-03-14 


54539.52 


-10.1 


3300 - 6030 


20.0 


— 




1 


2008-03-14 


54539.54 


-10.1 


5220 - 9240 


30.0 


N 


< 0.60 


1 


2008-03-19 


54544.55 


-5.1 


3140 - 9470 


0.5 


N 


< 0.06 


1 


2008-03-20 


54545.65 


-4.0 


3629 - 9437 


7.0 


N 


< 0.27 


1 








SN 2008bf 










2008-03-20 


54545.73 


-8.5 


3629 - 9437 


7.0 


A 


2.43 ±0.14 


1 


2008-03-24 


54549.61 


-4.8 


3805 - 10100 


7.0 


A 


1.48 ±0.05 


1 


2008-03-25 


54550.71 


-3.7 


3805 - 10214 


7.0 


A 


1.25 ±0.07 


1 


2008-03-26 


54551.71 


-2.7 


3805 - 10214 


7.0 


A 


0.78 ±0.07 


1 


2008-03-28 


54553.75 


-0.7 


3400 - 5300 


6.0 


N 




1 


2008-03-28 


54553.77 


-0.7 


4000 - 10200 


9.0 


N 


< 0.75 


1 
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TABLE 2 — Continued 



UT Date 


MJD 


Phase 


Wavelength 


Resol. 


C II 


piy(6300) 


Ref. 








Ran o-n ( A "1 
j_i,aiitie irt I 




V_/ldaa 








c?\ 
\ z > 


\.°) 


W 




V°) 


V 1 ) 










SN 2008bq 










2008-04-07 


54563.48 


0.3 


3800 - 9235 


8.0 


N 




i 








SN 2008fp 










2008-09-17 


54726.83 


-3.7 


3646 - 9460 


7.0 


A 


0.93 ±0.07 


l 


2008-09-21 


54730.85 


0.4 


3800 - 9235 


8.0 


N? 




i 








SN 2008gl 










2008-10-27 


54766.67 


-1.5 


3800 - 9235 


8.0 


N 




l 








SN 2008gp 










2008-11-02 


54772.73 


-6.0 


4025 - 10713 


6.0 


N 




l 








SN 2008hj 










2008-11-24 


54794.64 


-6.0 


3800 - 9235 


8.0 


N 


< 0.75 


7 


2008-11-25 


54795.57 


-5.0 


3800 - 9235 


8.0 


N 


< 0.63 


7 








SN 2008hv 










2008-12-10 


54810.81 


-6.6 


3300 - 6004 


27.0 






1 


2008-12-10 


54810.82 


-6.6 


5240 - 9210 


39.0 


\ 


0.93 ±0.25 


1 


2008-12-15 


54815.85 


-1.6 


3620 - 9429 


7.0 


F? 


< 0.18 


1 








SN 2008ia 










2008-12-10 


54810.77 


-2.8 


3300 - 6004 


2.7 


_ 




1 


2008-12-10 


54810.79 


-2.8 


5240 - 9210 


3.9 


F 




1 








SN 2009D 










2009-01-05 


54836.63 


-3.9 


3620 - 9428 


7.0 


N? 


< 0.24 


1 


2009-01-09 


54840.61 


-0.1 


5240 - 9210 


39.0 


N 




1 


2009-01-09 


54840.64 


-0.0 


3300 - 6004 


27.0 


_ 




1 








SN 2009F 










2009-01-05 


54836.66 


-5.4 


3620 - 9429 


7.0 


A 


3.58 ±0.15 


1 


2009-01-09 


54840.71 


-1.4 


5240 - 9210 


39.0 


N 




1 


2009-01-10 


54841.59 


-0.5 


3848 - 8017 


9.0 


N 




1 








SN 20091 










2009-01-17 


54848.56 


-5.0 


3629 - 9439 


7.0 


F 


1.35 ±0.26 


7 


2009-01-22 


54853.57 


-0.0 


3620 - 9428 


7.0 


N 




7 








SN 2009Y 










2009-02-08 


54870.83 


-5.5 


3716 - 9439 


7.0 


N 


< 0.24 


1 


2009-02-09 


54871.82 


-4.5 


3300 - 6004 


27.0 


_ 




1 


2009-02-09 


54871.83 


-4.5 


3716 - 9438 


7.0 


N 


< 0.18 


1 


2009-02-09 


54871.84 


-4.5 


5240 - 9210 


39.0 


N 


< 0.33 


1 


2009-02-09 


54871.84 


-4.5 


3694 - 6135 


2.0 


N 




1 


2009-02-11 


54873.87 


-2.5 


4032 - 10118 


5.0 


N 




1 








SN 2009aa 










2009-02-08 


54870.74 


-7.8 


3716 - 9439 


7.0 


F? 


< 0.18 


1 


2009-02-09 


54871.76 


-6.8 


3300 - 6004 


27.0 


_ 




1 


2009-02-09 


54871.78 


-6.8 


5240 - 9210 


39.0 


N 


< 0.81 


1 


2009-02-11 


54873.76 


-4.8 


4032 - 10119 


5.0 


N 




1 


2009-02-14 


54876.84 


-1.8 


4021 - 10119 


5.0 


N 




1 


2009-02-16 


54878.84 


0.1 


4021 - 10123 


5.0 


N 




1 








SN 2009ab 










2009-02-11 


54873.54 


-10.2 


4032 - 10123 


5.0 


N 




1 


2009-02-14 


54876.63 


-7.1 


4021 - 10120 


5.0 


N 




1 


2009-02-15 


54877.57 


-6.2 


4021 - 10125 


11.0 


N 




1 








SN 2009ad 










2009-02-15 


54877.59 


-8.3 


4021 - 10124 


11.0 


N 




1 


2009-02-16 


54878.63 


-7.3 


4021 - 10121 


6.0 


N 




1 


2009-02-23 


54885.59 


-0.5 


3800 - 9235 


8.0 


N 




1 








SN 2009ag 










2009-02-23 


54885.63 


-3.5 


3800 - 9235 


8.0 


N 


< 0.36 


1 


2009-02-24 


54886.70 


-2.4 


3800 - 9235 


8.0 


N 




1 


2009-02-26 


54888.68 


-0.5 


3800 - 9235 


8.0 


N 




1 








SN 2009cz 










2009-04-17 


54938.52 


-4.0 


3714 - 9435 


7.0 


N 


< 0.18 


7 








SN 2009dc 










2009-04-17 


54938.82 


-7.7 


3714 - 9433 


7.0 


SC 


19.66 ±0.10 


7 


2009-04-18 


54939.79 


-6.7 


3450 - 9658 


8.0 


sc 


18.86 ± 0.30 


7 


2009-04-22 


54943.78 


-2.8 


3360 - 9568 


8.0 


SC 


17.17 ±0.28 


7 


2009-04-23 


54944.80 


-1.8 


3361 - 9568 


8.0 


sc 


16.55 ± 0.29 


7 



NOTE. — Columns: (1) UT date of the observation; (2) Julian date of the observation (JD — 2400000); (3) Phase in days since S-band maximum 
light; (4) Wavelength range covered; (5) Spectral resolution in A as estimated from the FWHM of arc-lamp lines; (6) Type of 6300 A feature (see 
Sectio n |2. It ; (7) Pseudo-equivale nt w idth of absorption at £ s6300 (3-cr upper limits denoted wit h "<" sign); (8) Referenc es: 1 - IFolatelli et al 
| 20TT[) ; 2 - lAltavilla et alj 120071) ; 3 - IPastorello et al.| 120071) ; 4 - ITaubenberger et all 12008ft ; 5 - lYamanaka et all J2009a|) ; 6 - IStritzinger et al 
1 20101) ; 7 - This paper. 

f Spectrum corrected to match photometry fseelFolatclli ct al. |201lh . 
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TABLE 3 
Summary of SYNOW parameters 



Parameter 


SN 2004eo 


SN 2005el 


SN 2006ax 


SN 2009F 


Age [d] 


-6 


-7 


-10 


-5 


C II class 


N 


A 


F 


A 


v ph [km s _1 ] 


11,500 


12,300 


12,500 


12,000 


Tbb [K] 


12,500 


11,500 


12,500 


8,500 


r(C II) 






0.04 


0.1 


r(C II HV) 




0.05 






i> min (C II HV) [km s" 1 ] 




12,800 






r(C III) 






0.5 




no i) 


0.7 


0.5 


0.1 


1.5 


r(Na I) 


0.2 


0.1 


0.1 


0.5 


r(Mg II) 


1.0 


0.8 


0.7 


0.8 


r(Si II) 


4.3 


2.5 


3.0 


5.0 


r(Si III) 


0.(5 


0.5 


0.7 


0.8 


r S II) 


1.5 


1.0 


0.8 


0.5 


r(Ca II) 


14.0 


4.5 


7.0 


20.0 


r(Ca II HV) 


2.0 


1.0 


2.5 




iwCCall HV) [km s" 1 ] 


20,000 


20,000 


20,000 




r(Ti II) 


0.3 






1.0 


r(Fe II) 


0.3 


0.1 


0.1 


0.1 


r(Fe II HV) 


0.4 


0.2 


0.3 




iw(FeIl HV) [km s" 1 ] 


18,000 


19,000 


19,000 




T(Fe III) 


1.0 


0.5 


0.6 


0.5 



TABLE 4 

Fits of intrinsic (B — V) pseudo-color versus decline rate 



Sample 


a 


b 


CI 


NsNc 


(1) 


(2) 


(3) 


(4) 


(5) 


With carbon 


-0.023 ±0.028 


0.31 ± 0.15 


0.06 ±0.05 


6 


Without carbon 


0.020 ±0.018 


-0.01 ±0.07 


0.03 ± 0.04 


6 


All 


-0.002 ±0.018 


0.09 ±0.08 


0.05 ± 0.02 


12 


Folatelli et al. ( 2010b) 


-0.02 ±0.01 


0.12 ± 0.05 




14 



Note. — Fits of the type: (B — V)o = a + b[Am 15 (B) - 1.1] for Ami 5 (B) < 1.6 mag. 
Columns: (1) sample of SNe; (2) fit intercept; (3) fit slope; (4) intrinsic dispersion of fit (in magnitudes; see Section UJ; (5) 
number of SNe used in fit. 



